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We report on the diffusive-ballistic thermal conductance of multi-moded single-crystal silicon 
beams measured below 1 K. It is shown that the phonon mean-free-path 1 is a strong function of the 
surface roughness characteristics of the beams. This effect is enhanced in diffuse beams with 
lengths much larger than i, even when the surface is fairly smooth, 5-10nm rms, and the peak 
thermal wavelength is 0.6 pm. Resonant phonon scattering has been observed in beams with a pitted 
surface morphology and characteristic pit depth of 30 nm. Hence, if the surface roughness is not 
adequately controlled, the thennal conductance can vary significantly for diffuse beams fabricated 
across a wafer. In contrast, when the beam length is of order £, the conductance is dominated by 
ballistic transport and is effectively set by the beam cross-sectional area. We have demonstrated a 
uniformity of ±8% in fractional deviation for ballistic beams, and this deviation is largely set by the 
thermal conductance of diffuse beams that support the micro-electro-mechanical device and 
electrical leads. In addition, we have found no evidence for excess specific heat in single-crystal 
silicon membranes. This allows for the precise control of the device heat capacity with normal 
metal films. We discuss the results in the context of the design and fabrication of large-format 
arrays of far-infrared and millimeter wavelength cryogenic detectors. © 2014 AIP Publishing LLC. 
[http://dx.doi.org/10.1063/L4869737] 


I. INTRODUCTION 

The management of heat in dielectric beams and mem- 
branes is an important part of the design in micro-electro- 
mechanical systems ( MEMS). 1-3 While the physics of 
thermal transport in dielectric materials is well understood, 
the effect of fabrication processes on the surface physics is 
less clear and widely treated as a hidden variable in the eval- 
uation of the device performance. This uncertainty can pro- 
long the design-test cycle, where a parameter such as beam 
geometry is iterated until the target thermal conductance G is 
obtained. 

In our application, the MEMS is a cryogenic detector 
known as the Transition-Edge Sensor (TES). 6 Figure 1 
shows a TES designed for measurements of the polarization 
state of the Cosmic Microwave Background. 7 ' The TES is 
comprised of a superconducting MoAu bilayer deposited on 
a single-crystal silicon substrate. The silicon is etched around 
the bilayer and other electrical components to form a mem- 
brane that is supported by dielectric beams. The total thermal 
conductance G(T) of the beams is a critical parameter that 
determines the sensitivity (oc \/G), response time (cxl/G), 
and saturation power of the TES detector (oc J GdT) (satura- 
tion energy in the case of a TES micro-calorimeter). 6 The 
response time is also a function of the total heat capacity 
C(T) of the detector, which is the sum of the heat capacity of 
the silicon membrane and normal metal films. 
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At temperatures far below the Debye temperature of the 
lattice (0 S , = 645 K), anharmonic phonon collisions in the 
bulk of the dielectric are negligible, and the phonon-mean-free 
path £ is restricted by the cross sectional dimensions of the 
beam. In state-of-the-art far-infrared and millimeter TES 
detectors operating at sub-Kelvin temperatures, 4,5 ' 1 " the beam 
length L is much larger than £ , and the thermal conductance is 
often described by the Fourier law for phonon diffusion, 
G(T) = k(J) wt/L , where w and t are the width and thickness 
of the beam, and k(T) is the bulk thermal conductivity of the 
dielectric. k(T) is determined empirically with the fabrication 
and testing of numerous beams of varying width and length, 
and the aspect ratio h’/L is tuned to achieve the desired con- 
ductance. However, this approach rarely provides a complete 
view of the thermal physics, and higher order effects such as 
surface roughness of the beams have been identified as the 
likely cause of the reported non-uniformity in G across detec- 
tor arrays. 4,5 

A simple analysis can readily show the sensitivity of 
phonon propagation in long beams (L 3> £) to the details of 
the surface physics. For a random Gaussian surface, the 
roughness can be characterized by the ratio of the standard 
deviation a of the surface features to the peak thermal wave- 
length X th {T) of the incident phonon field. When a!X th -C 1, 
phonons are reflected specularly from the surface. 11 As the 
wavelength becomes much smaller than the scale size of the 
surface, a /l, /; 3> 1, phonons are scattered over all angles, or 
diffused. Assuming isotropic diffusion, it can be shown that 
for a rectangular beam, the phonon mean-free-path £ scales 
logarithmically with r = w/t (£ ~ 0.75 wr -1 In 2/- for r > I 
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FIG. 1. SEM image of a TES detector fabricated for a 40 GHz focal plane. 7 
The single-crystal silicon membrane volume is 5 x 350 x 450 /<m 3 . The ther- 
mal conductance from the membrane to the frame is effectively set by 
the ballistic-dominated thermal conductance of the short beam (magnified). 
The long silicon beams are5xl5x785 /mi 3 , and merely support TES bias 
and rf leads. The thermal conductance of the long silicon beams is 
diffusive-ballistic and sub-dominant to the conductance of the ballistic 
beam. The Pd layer deposited on the membrane sets the detector heat 
capacity. 


short beam, ballistic transport dominates and L/£ < 1. We 
emphasize that these limits bound the diffusive-ballistic con- 
ductance observed in practice. 

In Sec. II, we briefly describe the device fabrication and 
readout hardware employed for the low-temperature meas- 
urements of G and C. In Sec. Ill, we describe the theory of 
radiative phonon transport in diffusive-ballistic beams. In 
Sec. IV, we present results that illustrate the effect of fabri- 
cation on the diffusive-ballistic nature of the thermal con- 
ductance. In Sec. V, we describe the superior performance of 
ballistic beams for the control of G of TES detectors. 
Finally, we demonstrate that because of the low specific heat 
capacity of single-crystal silicon, the heat capacity of the 
detectors fabricated on this substrate can be effectively con- 
trolled using a normal metal film with high specific heat. 

II. METHODOLOGY 
A. Fabrication 


(Ref. 12)). Hence, in a beam where L> l, phonons experi- 
ence a large number of scattering events, N ~ LH, and the 
heat flux entering the beam at one end Pq(T) is reduced to 
P(T) = Pq(T) (1 - f) N at the exit port, where /is the diffuse 
fraction. The derivative dG/dN is a measure of the sensitivity 
of G to changes in N and proportional to the power-law func- 
tion ( 1 -ff. 

In practice, the statistics of the surface roughness can be 
heavily biased by the fabrication and is likely non-Gaussian. 
Dielectric beams can be roughened in various forms depend- 
ing on the compatibility and properties of the etchants with 
the dielectric material. Preliminary tests can shed light on 
the roughening steps, however, fabrication processes can be 
difficult to control in practice. Details such as base pressure 
and process chamber conditioning during dry-etch steps 
can vary between batches of devices, especially in shared 
facilities. As a result, the characteristics of the roughness 
of a single surface can change significantly over spatial 
scales shorter than the beam length, and between wafers. 
Furthermore, not all surfaces of a beam are rough to the 
same degree. These effects add to the unpredictable variabil- 
ity in N and therefore G in beams dominated by diffusive 
thermal transport. 

The work in this paper is focused on understanding the 
effects of fabrication on the thermal conductance in single- 
crystal silicon beams. The beam geometries explored span 
several orders of magnitude in L/wt. The results demonstrate 
that G in short (L S 1 0 /«n, L/£ < 1) uniform beams with 
ballistic-dominated phonon transport is insensitive to the 
fabrication conditions and can be realized with precision in 
devices on different wafers. This approach is suitable for 
TES detector arrays in ground-based and air-borne tele- 
scopes, where the detector noise equivalent power require- 
ment is greater than 1CT 18 W/a/Hz. 

Throughout this paper, the terms diffuse and ballistic are 
used to describe the scale size of the phonon transport in a 
silicon beam. In the presence of surface roughness, the ther- 
mal transport in a long beam with L/wt Jp- I is dominated by 
diffuse reflections off the beam walls, and L/i 1. For a 


The fabrication of the TES detectors reported in this 
paper is described in detail in Ref. 13, but for completeness, 
we provide a brief description here. Fabrication starts with a 
float-zone-refined p-type (100) boron doped silicon-on-insu- 
lator wafer. Heat flow is primarily in the [110] direction. The 
single-crystal silicon layer is 5 /mi thick, with a room tem- 
perature resistivity of 5kQ cm and a carrier density of 
5 x 10 11 cm -3 . A MoAu bilayer is deposited by a combina- 
tion of dc plasma sputter deposition and electron beam evap- 
oration. The Mo and Au are photolithographically patterned 
to define the superconducting sensor. The Au layer is etched 
by ion milling and the Mo layer is etched by reactive ion 
etching (RIE). Further metallic components are deposited on 
the device layer to enable the required operation of the pixel 
at rf frequencies, and to enable on-chip diagnostic capabil- 
ities such Johnson noise (JN) thermometry as needed. Au 
resistors are deposited by electron beam evaporation and pat- 
terned by lift-off. Following an in situ rf Ar plasma cleaning 
of the wafer, a Nb microstrip layer is sputter deposited, pho- 
tolithographically patterned, and etched in fluorine-based 
RIE chemistry. The silicon membrane is defined and reactive 
ion etched down to the buried oxide layer. The silicon wafer 
is then bonded to a pyrex handle wafer using a sacrificial 
wax. The backside of the silicon wafer is patterned by photo- 
lithography and the bulk silicon is etched by inductively 
coupled plasma deep reactive-ion etching (DRIE). The wax 
is then dissolved in solvent to release the devices. 

B. Experimental hardware 

Measurements were performed in a mechanically cooled 
system that reaches 3 K at the coldest stage. Further cooling 
of the detector stage down to 100 mK was achieved using an 
adiabatic demagnetization refrigerator. The devices were 
mounted and thermally anchored to a gold-plated copper 
package, and wirebonded to a Superconducting Quantum 
Interference Device (SQUID) cold stage amplifier for 
Johnson noise current measurements. Further amplification 
was achieved with a SQUID series array mounted at the 3 K 
stage. The temperature of the 100 mK cold stage was read 
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with a calibrated RuO, thermometer 14 with calibration trace- 
able to the National Institute of Standards and Technology. 
The precision of the temperature control was 32 /; K rms at 
150 mK. 

The room temperature electronics comprised of a cus- 
tom low-noise pre-amplifier, feedback electronics for flux- 
locked-loop operation, commercial low-noise high-gain 
amplifier and low-pass filter (SRS-650), and commercial 
16 Bit digital-to-analog electronics (NI PCI-6120) for digiti- 
zation of the amplified feedback signal. 

C. Johnson noise thermometry 

The method of JN thermometry is frequently utilized for 
the characterization of the thermal properties of thin-films 
below IK. 15 At audio frequencies (0.5-20 kHz) and sub- 
Kelvin temperatures, the current noise spectral density of a re- 
sistor R in series with the SQUID input coil takes the form 
Sj(f) = 4k B T e [1 + (f/f c ) 2 ]~ l /R , where T e is the electron tem- 
perature, f c = R/(2nL) is the characteristic frequency of the cir- 
cuit, and L is the total inductance of the input coil and stray 
inductance in the input circuit loop. This equation directly 
relates measurements of 5) to T e . For the test devices that are 
described in Ref. 13, the JN thermometer and heater are 
240 x 230 x 0.27 fi m 3 Au resistors deposited on top of interdi- 
gitated Nb leads. The Kapitza conductance at the Au-silicon 
interface and the electron-phonon conductance in the Au are 
approximately 10 nW/K each at 100 mK. 16,17 These conduc- 
tances are a factor of 50 greater than the highest silicon beam 
conductance measured in our devices. We, therefore, assume 
the electron system in the thermometer is in thermal equilib- 
rium with the phonons in the silicon membrane. The maxi- 
mum expected temperature gradient across the membranes is 
less 0.5 mK across the 0.1-1 K temperature range. 

The random error in a temperature measurement from 
the JN current can be estimated from the radiometer equation 
o T = T /y/B z m , 15,18 where 5 = 5 kHz is the bandwidth over 
which the current noise power spectral density was measured 
and z m = 102 s is the measurement time. Hence, the frac- 
tional error throughout the temperature range of our meas- 
urements is less than 0.5%. Measurements of the JN current 
as a function of bath temperature were used to fully calibrate 
the cryogenic and warm readout amplifiers. Typical readout 
noise of 15pA/\/Hz referenced to the input circuit was 
achieved. The JN signal was 30 — 100pA/\/Hz in the range 
of 0.1-1 K. 

To measure G of the silicon beams, a small signal dc 
heater power was applied to the membrane and the increase 
in temperature was recorded. The relationship between 
heater power and membrane temperature at a given bath 
temperature T B can be parameterized by the power law func- 
tion P = K(T" — Tg), where P is the heater power, T is the 
membrane temperature, K is a constant that is proportional 
to the number of phonon modes, and n is the temperature 
exponent that is indicative of the temperature dependence of 
the phonon mean-free-path. The first order derivative yields 
the thermal conductance G(T ) = nKT"~ l . In our measure- 
ments, the change in temperature was small, ~10%, and the 
approximation G(T) = AP(T) /AT was valid. 


D. TES current-voltage thermometry 

We have also utilized TES current-voltage curves to 
determine the thermal conductance of detectors in which 
normal metal films for JN thermometry were omitted. When 
a TES is voltage-biased on the normal-to-superconducting 
phase transition, the silicon membrane temperature is 
fixed at the critical temperature of the MoAu bilayer, T c 
= 150 mK. The thermal conductance is then determined 
from a fit to the Joule power dissipated in the TES as a func- 
tion of bath temperature. In this case, P = Pj = KiT'/ — T'/ ) 
and G(T b ) =nKT B ~ l , where the parameters K and n are 
equivalent to those derived from JN thermometry in the 
same temperature range. 

III. MODELING HEAT TRANSPORT 

In building a model of heat flow in a silicon beam, it is 
important to first identify the thermal baths that define the 
boundary conditions in the model. A thermal bath is ideally 
a blackbody that provides a matched termination for phonons 
at all wavelengths. In a detector such as that shown in Fig. 1, 
each silicon beam is connected to two thermal baths, the 
frame and the membrane. The electron system in the Pd film 
is strongly coupled to the phonon system in the membrane. 
At 150 mK, the electron-phonon conductance in the Pd is 
30 nW/K, 1 ’ and the Kaptiza conductance between the Pd 
and the silicon membrane is ~ 1 30 nW/K. k> For comparison, 
the highest beam conductance measured is 240 pW/K at 
150 mK. The silicon frame is effectively a blackbody since 
its lateral dimensions are much larger than the peak thermal 
wavelength. At 100 mK, the peak thermal wavelength is X th 
~ Ahc/k B T = 0.6 /tm in silicon, where h and k B have the usual 
definitions and c = 5070 m s 1 is the speed of sound for 
transverse waves. A phonon scattered from a beam into the 
silicon frame has a low probability of coherently backscatter- 
ing into the beam. Given the 5 nm thickness of the silicon 
layer, and the thermal wavelengths above 100 mK (A f * w 0.6 
(0.1/7) /fin), the membrane and the frame can be modeled as 
three-dimensional multi-moded thermal baths, where the 
number of modes is given by the Planck distribution and is 
much larger than unity. 

The beams behave as three-dimensional multi-moded 
acoustic waveguides since the widths (w > 13 /un) and 5 /an 
thickness are much larger than X tll at 0. 1-1 K. Phonons emit- 
ted by the baths couple radiatively to the silicon beams, 
propagate unimpeded in the bulk of the beams, and scatter at 
the beam surfaces. When the surfaces of the beams are per- 
fectly smooth, phonons coupled to the beams are scattered 
elastically without any loss of coherence in the beams. The 
conductance in this case is determined by the beam areas 
and is referred to as ballistic. In the case of rough surfaces, 
phonons scatter diffusively with a large probability of back- 
scatter. The scattering is elastic but the phonon wavevectors 
are randomized at the surface. If the surface features form 
cavities much larger than A,*, the surface can be thought of 
as a blackbody. The conductance is thus boundary-limited 
and was first described theoretically by Casimir. 211 An impor- 
tant assumption in Casimir’s formalism is the isotropic 
(Lambertian) re-emission of phonons by the beam surfaces. 
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In practice, the scattered pattern of phonons from a 
rough surface is strongly dependent on the surface character- 
istics. When the roughness can be described statistically, the 
reflection pattern is averaged over the ensemble of the sur- 
face realizations. For a surface with a Gaussian correlation 
function, the roughness can be characterized by the correla- 
tion length £ and standard deviation a of the surface features 
above a mean height. Figure 2 illustrates the scattered pattern 
as a function of er/A. 1 1 For a nearly smooth surface, <r/A <§; 1 , 
phonons are largely reflected in the specular direction. When 
rr/z 1 , phonons incident on the surface are diffused in all 
directions, however, the flux is not Lambertian. Nonetheless, 
in numerous theoretical studies of boundary-limited 
scattering, ’ the reflected pattern is decomposed into a 
diffuse Lambertian fraction / and a specular fraction 1 -/, as 
illustrated in Fig. 2. The diffuse reflection is modeled as 
inelastic scattering that redistributes the phonon wavevector, 
although strictly speaking all reflections are elastic at the sur- 
face facets that define the roughness. The surface can then be 
thought of as a blackbody with emittance /. Hence, a com- 
pletely lossless but rough surface can be modeled as a lossy 
surface. ’ " This approach is numerically less intensive 
than calculations of the reflection pattern based on the exact 
surface profile, or the statistics of the surface. 1 1 

Our model of heat flow in a silicon beam has the follow- 
ing assumptions: (i) The ends of the beam constitute thermal 
baths with constant temperature and unity diffuse fraction. 
The baths continue to distances much larger than the beam 
length as shown in Fig. 3(a). This construct is numerically 
convenient from the perspective of the heat transport model 
presented, and is similar to that adopted in Ref. 22. (ii) The 
surfaces of the beam reflect phonons with a diffuse fraction/ 
that is isotropic, and a specular fraction given by 1 -/. (iii) 
The gradient of the temperature profile across a beam is line- 
arly dependent on/ and constant as a function of length. This 
assumption is valid in multi-moded beams, where the num- 
ber of modes is much larger than unity and the diffuse pho- 
non reflection can be treated as diffuse emission. _n Models 
that iteratively solve for the temperature profile in multi- 
moded beams show negligible deviations from a constant 



FIG. 2. Beam patterns of phonons reflecting off a fairly smooth {a IX <C 1) 
and a very rough surface (a/X ^>1). The diffuse fraction of the realistic 
beam patterns is generally approximated as isotropic (Lambertian), which is 
relatively simple to integrate into models of thermal transport. 


(a) 


L 


-> 




FIG. 3. (a) Beam geometry and linear temperature profile assumed in the 
model of heat transport. The temperature profile is linear across the beam. 
The temperature is constant at the thermal baths, which are represented by 
diffuse surfaces (f= 1) that extend to distances much larger than L. (b) Each 
surface element of the beam radiates phonons into all solid angles. The radi- 
ated flux from the infinitesimal area dAj that reaches the shaded beam port is 
reduced depending on the diffuse fraction of the beam surfaces, (c) Using 
the method of images, the path of specularly scattered phonons can be re- 
imaged to one vector. The components of the vector map the number of scat- 
tering events at each surface of the beam. 

temperature gradient. 2 ' When the phonon mean-free-path is 
much smaller than the beam length, or / is close to unity, 
phonon diffusion is the dominant term of the total heat flux, 
and the constant temperature gradient is an excellent approx- 
imation. The approximation also holds when / approaches 
zero, or the mean-free-path is comparable to the beam 
length, since the temperature gradient is linearly dependent 
on / (the mathematical description of this detail is in the 
Appendix), (iv) The phonon modes that propagate are two 
orthogonal transverse modes with a sound speed of 
v,= 5070 m/s, and one dilatational mode with a sound speed 
of v t = 13 600 m/s. Each mode has a thermal spectrum, (v) 
Phonons propagate unimpeded in the crystalline bulk of the 
silicon beams, and only scatter at the beam surfaces, (vi) The 
diffuse flux is re-emitted according to Lambert’s law of ther- 
mal radiation, (vii) All modes scatter diffusively with the 
same probability/, and there is no spectral dependence for/ 
It is possible to map /(A,/,) by fitting measurements of 
G(T) and calculating the peak thermal wavelength as a func- 
tion of Tr For a Gaussian surface," when er/A <C 1, the 
specular power reflection coefficient is proportional to A -2 
and the diffuse coefficient to A~ 4 ; when <r/A S> 1, the reflec- 
tion coefficient is independent of wavelength. Hence, the 
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thermal spectrum can be scattered very differently as temper- 
ature changes, and in recognizing the limitations of our 
model, we do not extrapolate the results to f(X tk ). 

Figure 3(b) illustrates the radiative model. The beam is 
divided into infinitesimal areas clA, = dy t dx, that radiate heat 
into all solid angles. The phonons are scattered several times 
before reaching the beam exit port and thermalizing in a 
bath. During each scattering event, the phonon flux is 
reduced by a factor f k that is the diffuse fraction of the sur- 
face k. The nature of the problem is then to simultaneously 
calculate the flux that is viewed by dAj from dA, and count 
the number of scattering events from each surface of the 
beam. For the latter, the complex path of the phonons can be 
unfolded, or re-imaged, to a single vector that reaches an 
equivalent beam exit port as shown in Fig. 3(c). The compo- 
nents of the vector then equate to the number of scattering 
events from each surface. 

The details of the radiative model are described in the 
Appendix. We define an equivalent view factor for a rectan- 
gular waveguide 

4 r 2 

F{xi) = ^2 F^(f m ;Xj) + y ^F ± (s,f m -,x i ) +F D (f m ;Xi) , 

m= 1 [_ 5=1 

(i) 


where F“(f m ',xi) and F ± (s,f m \x j ) are equivalent view factors 
between parallel and perpendicular surfaces of the beam, 
respectively, and F D (f m \xi) is the contribution of the direct 
(unscattered) flux from the surfaces of the beam to the cross 
sectional area at a = 0 (reference plane). f m is the diffuse 
fraction of the surface emitting the flux, and the index m runs 
over all surfaces of the beam. The index s runs over the 
surfaces perpendicular to m. Included in the terms F^ (f m :xj) 
and F ± (s,f m -,Xi) is a factor n!=i (1 ~fk) Nt , where N k is the 
total number of scattering events from each surface as 
described in Fig. 3(c). For the baths, F B (xj) is obtained by 
setting/,,, = 1 in Eq. (1). The total heat flux through the beam 
isP = P\ x= q, 


p — q 


T A [xi)F[x,) 

J-L/2 


f L/2 


J-L/2 

r-L 


T*(xi)F(xi ) 


+/ 


F B (xi) - T\ 


F B {xi) 


dx; 


(2) 


where q = ti 2 k B / (I20ti 3 v 2 ) , and v = y 2/v 2 + l/vf is an av- 
erage group velocity. Given the linear temperature profile 
T{xi) = Xj ( dT / dxj ) + T c , the approximation T k — T* « 4 
(. dT/dxj)LT 3 , where T is the average of the bath tempera- 
tures, and the relation F( x,-) = F(— x,), Eq. (2) reduces to 


P = 4q T 3 — ( 2 

dxi \ 

dT 


■L/2 


X{F(xi) dx, + L 


F B {xj) dxj , 


= 4qT i - r (2I d + LI s ), 
dx. 


( 3 ) 


where I d is the net diffuse power from the beam surfaces, 
and / is the net specular power that is exchanged between 


the baths. The conductance of the beam as a function of tem- 
perature when scaled by the maximum (ballistic) conduct- 
ance G/j = 4qT 3 A is 


8 


G _ ( P A\ 1 
Gb~ \ A{dT/dxi )l)Gb' 




( 4 ) 


In the short beam limit, / d — > 0 and I s — > A. For non-zero dif- 
fuse fractions and as the beam length becomes much larger 
than the cross-sectional dimensions, I s —* 0, and using 
k = Cyvlj 3, where Cy = \6qT 3 /v is the volume specific 
heat of the material, the mean free path takes the form 
i = 3I d /(2A). 

When the diffuse fractions of the surfaces of the rectangu- 
lar beam are different, all terms in the summation over the 
index m in Eq. (1) must be evaluated. The numerical evalua- 
tion of the problem is reduced by a factor of two if there is 
symmetry in the diffuse fraction of the surfaces, /r op = /Bottom 
and /Left = ./kight- For the silicon beams discussed in this paper, 
the top surface and sidewalls are diffuse, and the bottom sur- 
face is specular. It is important to note that the diffuse frac- 
tions of the top and sidewall surfaces are uncorrelated. 

The terms F"(f m ;xi) and F- L (i,/ m ;x,) in Eq. (1) are eval- 
uated using Monte-Carlo integration with 10 5 uniformly dis- 
tributed random points averaged 10 2 times. F° (/„, ; x, ) is 
evaluated with an adaptive quadrature routine. Figure 4 illus- 
trates the accuracy of the model compared to an analytic 
expression for the mean-free-path i in a rectangular beam. 1 " 
At the diffuse limit ( L/£ 'J> 1), the conductance scales as 1/L, 
and at the ballistic limit, the conductance is independent of 
length. The expression g = H(L + £) is remarkably accurate 
over the range of beam length explored in Fig. 4. This fact 
highlights the importance of the phonon mean-free-path, that 
a single parameter can capture a significant fraction of the 
heat transport in the ballistic-diffusive limit. 



FIG. 4. Simulations of a beam with fixed cross-section compared to the ana- 
lytic expression in Ref. 12. The sidewall diffuse fraction is unity. The stand- 
ard error from the Monte-Carlo integration is shown for each simulation. 
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IV. DIFFUSIVE-BALLISTIC CONDUCTANCE 

Figure 5 shows the thermal conductance of silicon 
beams as a function of beam area-to-length ratio. The cause 
of the variability in conductance (a factor of ~5 for the same 
beam geometry) is the roughening of the silicon beam surfa- 
ces during fabrication. The top surface of the beams is pri- 
marily roughened during a dry etch step that patterns a Nb 
layer. The plasma chemistry, 75% CF4, 15% 0 2 , etches sili- 
con at a higher etch rate than Nb. The exact processes that 
roughen the silicon are unclear, but the likely candidates are 
micro-masking due to non-volatile species in the plasma, 
and slower etch of native and in situ NbO, in the fluorine 
chemistry. The use of an etch stop such as lOOnm thick 
A10 v was successfully employed in early devices; however, 
these detectors exhibited very large uncontrolled heat 
capacity due to the amorphous nature of the etch stop. 24 Test 
devices were also fabricated with a lift-off process that left 
the beam surfaces smooth. The conductance of these beams 
was higher, suggesting boundary-limited phonon propaga- 
tion as the cause of the variability in conductance of the 
beams with rough surfaces. 

We have explored the relationship between thermal con- 
ductance and surface roughness with an extensive dataset of 
atomic force microscopy (AFM) images. Non-destructive 
AFM of a silicon beam surface a millimeter in length is not 
practical. Instead, the silicon device layer was imaged on 
the handle wafer, in close proximity to each device. For a sil- 
icon surface with a Gaussian roughness, a ~ 5-20 nm and 
C ~ 100-200 nm depending on the processing conditions. 
However, we have found that a large fraction of the rough- 
ened surfaces did not satisfy the mathematical assumptions 
required to apply the theory of phonon scattering by nearly 
smooth Gaussian surfaces, and the measured thermal 
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FIG. 5. Conductance of single-crystal silicon beams as a function of beam 
area to length ratio. The conductance can vary significantly depending on 
the fabrication conditions. Devices that were dry etched during fabrication 
had roughened silicon surfaces. The dashed line is the fit G = kA/L (exclud- 
ing the lift-off points), and is the simplest interpretation of the data that is 
intended to merely guide the eye. Details of this data are summarized in 
Table I. 


conductance differed from expectations of such a theory, 
even when a was much less than the typical phonon wave- 
length at the temperatures explored (0. 1-1 K). We observed 
power spectral densities of surface height that were approxi- 
mately power laws over the range of observable wavevector 
values. In this case, we found an empirical relationship 
between the power-spectral-density S(k) of the surface profile 
and the measured magnitude and temperature dependence of 
the phonon mean-free -path in our silicon beams, 
t = nt / (4k 2 \/S(k) / (27r) ) , where t is the silicon beam thick- 
ness and k= 0.010 T/(nm K) is the phonon wavevector. As 
expected from Beckmann’s theoretical treatment of scattering 
from non-Gaussian surfaces, 24 the observed relationship 
involves the correlation function of the local tilt angle of the 
surface facets, 0 rms (k) = k 2 y/S(h)/(2n). 

It is challenging to precisely predict the phonon mean- 
free-path from information contained in the AFM data due 
to practical sampling considerations. This occurs for sev- 
eral reasons: The statistics of the surface roughness are 
observed to vary significantly across the wafer and are 
sparsely sampled in 8 x 8 /rm 2 patches set by the resolution 
and stability of the AFM. In addition, for a multi-moded 
beam with thickness of order the width, details such as sur- 
face roughness of the beam sidewalls are significant, as is 
shown in Sec. V. Non-destructive AFM imaging of the 
sidewall geometry is extremely difficult to perform, and 
such a task becomes impractical when considering the 
number of silicon beams that constitutes a detector imaging 
array. However, we have found that AFM images can 
provide qualitative insight to the thermal conductance in 
beams where the surface roughness relative to the 
thermal wavelength has a significant amplitude and is 
non-Gaussian. 

Figure 6(a) shows the conductance of smooth (lift-off) 
and roughened (dry etch) silicon beams when compared to the 
boundary- limited conductance in each case. 12 The phonon 
propagation transitions from diffuse to ballistic at 300 mK, 
where the peak thermal wavelength is /. th ~ 200 nm. The min- 
imum in conductance suggests the transition is resonant for 
the rough beam. 2 ' The inset in Fig. 6(a) shows the inverse 
phonon mean free path 1 . The scale height of the scatterers 
C can be estimated from the temperature at which £~ 1 peaks. 
Using the Rayleigh criterion, £ « 8 sin y, 1 1 where /,/, is 

the phonon wavelength at the peak temperature, and y is the 
phonon angle of incidence to the surface. Assuming an aver- 
age over all angles, (siny) =2/71, the Rayleigh criterion sug- 
gests C ~ 40 nm for the rough beam and ( S 6nm for the 
smooth beams. These estimates are consistent with the rough- 
ness measured from the AFM images shown in Fig. 6(b). The 
rough (dry etch) surface is heavily pitted, with an average pit 
depth of 30 nm and average size of 250 nm that is comparable 
to A th at 300 mK. The scale size is much smaller for the 
smooth (lift-off) surfaces, a ~ 3 nm, which leads to the onset 
of resonance at a much higher temperature. 

In addition to phonon scattering by rough surfaces of a 
beam, organic residue such as photoresist on the surfaces can 
enhance diffuse scattering. Thin layers of metals and non- 
metals hundreds of Angstroms thick have been shown to 
induce diffuse transport in silicon. 26 In our devices, the 
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FIG. 6. (a) Conductance of silicon beams scaled by the boundary-limited 
conductance in a rectangular beam (dashed line). 12 Solid lines are spline fits 
to guide the eye. The labels refer to the top surface of the beams, which is 
roughened by a dry etch step during fabrication (squares), but is left smooth 
when the dry etch step is replaced with a lift-off process (diamonds, trian- 
gles). The inset shows the inverse phonon mean-free-path as a function of 
temperature. The peak in i~ x suggests resonant phonon scattering from the 
rough surface, (b) AFM images of 6.4 x 6.4 /mi 2 patches representative of 
the top surface of the beams described in (a). The scale size estimated from 
the peaks in r is consistent with the roughness scales measured from these 
images. The rough surface is heavily pitted, with an average pit depth of 
30 nm. The large distribution in pit size and shape is the likely cause of 
the observed broad resonance. Since the scale size is much smaller for 
the smooth surfaces ( o ~ 3nm), the onset of the peaks occurs at a much 
higher temperature (shorter thermal wavelengths). The striping in the AFM 
image of the smooth surface is due to residue and accumulation of particu- 
lates on the AFM tip. Beam dimensions are L={0.4, 1.5, 3} mm, and 
w = { 24, 20, 25 } //m for devices labelled with squares, diamonds, and trian- 
gles, respectively. 

reduction and variability in G in a few cases were at least 
qualitatively attributed to residue accumulation on beams of 
order a millimeter in length. Care must be taken when devi- 
ces are cleaned since the cleaning process may introduce fur- 
ther residue, or roughness (see example of silicon dipped in 
HF in Ref. 26). 

Table I summarizes the silicon beam geometries and 
conductances shown in Fig. 5. The exponent ft = n - 1 is 
inversely proportional to G, indicating a transition from dif- 
fuse to ballistic transport. Low G devices have boundary- 
limited transport, and /) ~ 3. High G devices are approaching 
ballistic transport, and /? < 3 in the transition region. 


TABLE I. Conductance fit parameters a = nK and /> = n - I as a function of 
beam geometry, G = 77’ . The entries marked with an asterisk indicate 
beams with smooth top surfaces (lift-off devices described in the text). 


Length 
in m) 

Width 

(/mi) 

ot (pW/K^ +1 ) 

P 

G (150 mK) 
(pW/K) 

5000 

35 

1054 ± 55 

3.13 ± 0.06 

2.78 

5000 

45 

1449 ± 50 

2.96 ± 0.04 

5.28 

3000 

25 

499 ± 36 

2.92 ± 0.08 

1.96 

3000 

25 

1241 ± 118 

3.27 ± 0.09 

2.51 

3000 

25 

1210 ± 103 

2.19 ± 0.05 

18.99* 

1500 

15 

625 ± 44 

2.92 ± 0.09 

2.46 

1500 

15 

583 ± 31 

3.21 ± 0.07 

1.32 

1500 

20 

2000 ± 81 

2.38 ± 0.03 

21.88* 

1000 

25 

2406 ± 37 

2.89 ± 0.03 

10.00 

553 

20 

2521 ± 67 

3.11 ± 0.03 

6.90 

544 

20 

1922 ± 132 

2.80 ± 0.08 

9.48 

400 

24 

4292 ± 386 

2.45 ± 0.06 

41.10 


V. CONTROL OF G IN SILICON-BASED DETECTORS 

In the design of prototype TES detectors, the choice of 
the beam length and width was guided by the data in Fig. 5. 
Several strategies were adopted in order to control the con- 
ductance and produce uniform detectors for the focal plane 
array described in Refs. 7 and 8. To minimize the top surface 
area of the silicon beams roughened by the dry etch, the 
width of the Nb leads was extended to cover a larger fraction 
of this surface, and the beam widths were defined with a 
combination of DRIE using the standard Bosch process 27 
followed by RIE to remove the silicon down to a Nb ground 
plane. The two processes are separated to reduce potential 
metallic contamination in the DRIE chamber, while mini- 
mizing lateral silicon undercut from the isotropic RIE. The 
bottom surface of the beams in contact with the Nb ground 
plane is smooth, and scanning electron microscopy images 
of the sidewalls showed no significant features from the RIE 
etching of the beams employed in test structures, Fig. 7(a). 
The beam sidewalls were predicted to be nearly specular, 
with conductances comparable to the lift-off devices shown 
in Fig. 5. The top surface of the beams was expected to be 
slightly rough, /~ 0.2, since the area exposed during the dry 
etch step was reduced by a factor of two. Four beams 

(a) RIE Si etch (b) x3 DRIE + RIE Si etch 


FIG. 7. SEM images of silicon beams. The Nb lead is 200 nm thick, (a) The 
isotropic RIE chemistry leaves smooth sidewalls with a slight taper. In 
Fig. 5, devices marked with crosses, squares, diamonds, and circles have 
similar sidewall profiles, (b) To minimize the top surface area roughened by 
the dry etch of the Nb lead layer, DRIE was used to define the beam width. 
The ridges at the top of the beam are scalloping from the DRIE cycles. 
Possible sources for the bright features on the sidewall are leftover passiva- 
tion layer and charged silicon features. 
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785 /mi in length and 13 /im in width were integrated into the 
TES design. 

Post fabrication, the beams had well-defined top surfa- 
ces, but the sidewalls were roughened by the DRIE steps, 
Fig. 7(b). The total conductance was a factor of 5 below the 
predicted target range. In principle, the scale length and 
height of the sidewall features that include the DRIE scallop- 
ing as well as passivation residue present a very diffuse sur- 
face for phonons above 100 mK. Figure 8 compares the total 
conductance of the four beams to simulations based on the 
diffusive-ballistic model described in Sec. III. The diffuse 
fractions of the sidewall and top surface of the beams vary 
across a wafer and between wafers. 

The inferred diffuse fraction of the top surface of the 
beams is f < 0.4. A possible source for the moderate rough- 
ness is the 200 nm thick Nb lead sputter-deposited on each 
beam. The lead and beam are generally considered as 
decoupled systems, and the Nb conductance is estimated to 
be far below that of the silicon beam. 1 However, it is the dif- 
fusion process of the phonons in the silicon by the Nb layer 
that is of importance and generally unknown. It has been 
shown that even in the presence of large acoustic mismatch, 
phonons in silicon can be diffused at surfaces with metallic 
layers tens to thousands of Angstroms thick (see Fig. 8 in 
Ref. 26 for examples of phonon diffusion caused by A1 and 
Ti films below 1 K). 

It is evident from the results in Fig. 8 that a solution to 
the problem of conductance control is a design that is least 
sensitive to the effect of fabrication conditions such as sur- 
face roughening and residue. Our approach was to reduce the 
length L of a beam, which increases the phonon mean-free- 
path t by reducing the interaction probability of the phonons 
with the rough sidewalls of the beam. The conductance is 
less sensitive to £ in the limit Ll£ < 1. Figure 9 shows the 
simulated and measured total conductances when a ballistic 
beam 10 /un long and 1 3 /mi wide is integrated into the 




FIG. 8. (Left) Conductance of the four 785 /m 1 long beams shown in Fig. 1. 
The measured conductance was a factor of ~5 below the design conductance 
(hatched box). The cause of the variability is boundary-limited scattering. 
(Right) Simulations of the conductance as a function of beam wall diffuse 
fractions. The inferred diffuse fraction of the samples are moderately rough 
top surfaces, / Top < 0.4, and very rough sidewalls, /sidewalls > 0.6. 
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FIG. 9. (Left) Simulations of the total conductance of the beams including 
the short ballistic beam shown in Fig. 1 . The shaded and hatched regions are 
defined by the measured range of conductance of the long diffuse beams 
shown in Fig. 8. (Right) The measured conductances are within the predicted 
range of the simulations (light grey). 


detector as shown in Fig. 1 . Given our knowledge of the sur- 
face roughness after the dry etch step (Fig. 8), a diffuse frac- 
tion / < 0.4 was expected and subsequently obtained in 
practice. 

The total conductance, which is largely set by the ballis- 
tic beam, was predicted prior to fabrication. The diffuse 
beams merely support microstrip and bias leads necessary 
for the operation of the detector. However, they contribute 
the largest fraction to the variance in the total conductance 
across a wafer and between wafers. The magnitude and var- 
iance in G of these beams were estimated from the data in 
Fig. 8. The fractional standard deviation in total conductance 
is ±8% and within the requirement of the target instrument. 7 

The range of conductances that can be accessed with a 
ballistic beam is suitable for low-noise bolometric detectors 
operating at ground-based or air-borne telescopes where the 
signal power is of order 1-3 pW/K. In the ideal design, a sin- 
gle ballistic beam sets the conductance of the detector and 
carries all signal and bias leads. This approach is very well 
suited to the fabrication of kilo-pixel arrays, especially when 
the signal is coupled to the detector with free-space absorb- 
ers. For space applications where the signal power is less 
than 1 pW, the ballistic beam can be a phononic structure 
with phonon bandgaps that can be engineered to reduce the 
conductance. 22s 


VI. CONTROL OF CIN SILICON-BASED DETECTORS 

The heat capacity of the dielectric also plays an impor- 
tant role in the performance of thermal detectors. For a bo- 
lometer and calorimeter, C affects the detector time constant, 
and in some cases anomalous heat capacity associated with 
the dielectric has been reported to increase thermodynamic 
noise of these detectors. 29 

The heat capacity of the silicon membranes was deter- 
mined from C = Gi lh , where the thermal time constant x th 
was measured from the decay in temperature of a silicon 
membrane to a pulse in heater power. Compared to the 
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Debye specific heat, 9 Cy ebye /Cy eas = 0.2 at 200 niK and 0.7 
at 1 K. The systematic difference between the measurements 
and the Debye approximation is due to the large contribution 
to the total heat capacity of the devices from the normal metal 
Au heater/thermometer and Pd-Au termination on the silicon 
membranes. Nonetheless, the measured specific heat is 
extremely small. For the detectors shown in Fig. 1, the esti- 
mated silicon heat capacity is ~2 fl/K at 150 mK, which is ad- 
vantageous since the overall detector heat capacity can be 
controlled with a normal metal that can be reliably deposited 
on the membrane. We have successfully used Pd to achieve a 
heat capacity of 2 pJ/K. 7 Measurements indicate that the spe- 
cific heat of Pd below 1 K is 5.8 mJ mol 1 K 2 . which is lower 
than the room temperature value of 9.4 mol 1 K 2 . We have 
successfully deposited 400 nm Pd layers with a lift-off process. 

It is important to emphasize that there is no evidence for 
anomalous specific heat in the single-crystal silicon device 
layers. This is in contrast to the large excess specific heat of 
amorphous SiN v that is widely used in the fabrication of ther- 
mal detectors. 29-3 1 


the NASA Postdoctoral Program at the Goddard Space Flight 
Center duting the initial stages of this work. We thank Samelys 
Rodriguez for wirebonding and cryo-cable support. We thank 
our collaborators A. Ali, J. Appel, and T. A. Marriage for their 
contribution to the data acquisition. We thank the National 
Institute of Standards and Technology for providing the 
SQUID series array and time division multiplexer chips. 

APPENDIX: RADIATIVE PHONON TRANSPORT 

We start from the usual definition of the infinitesimal 
view factor dF c [ A ,-dA l ,’' 2 

cos 6; cos 9 j 

dAj dFjA^ = -=T2 — dAj dAj , ( A 1 ) 

n\R\~ 

where R is the vector between the infinitesimal surfaces dA, 
and dAj , and 0 V are the angles between the normal vectors to 
dAij and R (see Fig. 3). We define an equivalent view factor 
for convenience, 


VII. CONCLUSIONS 


f 7 dAj— dAj — dAj dF dAi- dAj , (A2) 


The thermal conductance of single-crystal silicon devi- 
ces can be precisely controlled using a short beam with 
ballistic-dominated phonon transport. The conductance of 
the short beam is largely set by its area and is insensitive to 
the details of the fabrication conditions that can vary over 
time. This approach is superior to the conventional use of 
long beams for the control of thermal conductance. Phonon 
transport in long beams is diffusive-ballistic and very sensi- 
tive to the detailed surface physics of the beams. The con- 
ductance is thus boundary-limited, and has been observed in 
beams with modest surface roughness (5-10 nm rms in sur- 
face height). Resonant phonon scattering has also been 
observed at 300 mK in beams with highly pitted surfaces 
with pit depths of order 30 nm. Boundary-limited scattering 
has contributed to a large variance, up to a factor of 5, in the 
conductance of devices fabricated and tested during this 
work. In contrast, when a beam with ballistic-dominated 
thermal transport is integrated into the design, the uniformity 
in conductance is reduced to a fractional deviation of ±8% 
in devices fabricated across two wafers, and this variability 
is largely set by the conductance of diffuse beams that sup- 
port electrical and microwave leads to the device. In addi- 
tion, we have found no evidence of excess specific heat in 
single-crystal silicon membranes. Hence, the total heat 
capacity of the detector can be effectively controlled with a 
normal metal film. For the TESs described in this work, the 
heat capacity was determined by a Pd layer of 400 nm thick. 
Hence, the strategies outlined for the precision control of 
thermal conductance and heat capacity are well suited to the 
fabrication of uniform large-format arrays of TESs with sen- 
sitivities approaching 10 -18 W/Hz 1,/2 . 


and note that the heat exchanged between the infinitesimal 
surfaces is 

dPdAi-dAj(x) = q [Tf(x) - Tf(x)} F dAi- dAj, (A3) 

where Tjj(x) are the temperatures of the surface elements. 
The total heat flux is the integral of dPdA t -dAj over the beam 
length. 


P = 


dP dA,— dAj (-v 5 
. S 


(A4) 


where S represents the surfaces of the waveguide, which in 
our case is a rectangular beam. 

The equivalent view factor between any two perpendic- 
ular surface elements of the beam is 


dAj— dA; 


yizj 


n[(xj-Xj) 2 +zf+yf ] 2 




,N k 


„ dAjdA 


7> 


=J r - L (m-,Xi,y i ,y j ,z j )dx i dy t dyj dzj , 


(A5) 


where the indices m and k sum over all surfaces of the beam. 
For parallel surface elements, 


r 


dAj— dAj 


n[{ Xi - Xjf + {y/ - yj) 2 + t 2 } 2 


/.IK 1 -/*) 


N t 


dAi dAj, 


= T^{f m -,Xj,yi,yj,Zj) dxidytdxj dyj. (A6) 
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N k is the number of specular reflections experienced by pho- 
nons reflecting off the element dA, first, and then all other 
surfaces until a reference plane (e.g., x=0 in Fig. 3(a)) is 
reached. N k can be determined from the geometry of the 
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specular scattering as described in Fig. 3(c), and is the sum 
of scattering events from each surface. For example, if the 
radiating element is on the top surface, N k can be found for 
phonons reflecting off the top and bottom surfaces, and 
reflecting off the sidewalls only. In this way, a different dif- 
fuse fraction can be assigned to each surface. It is also useful 
to note that the linear dependence of the equivalent view fac- 
tors on f m can be translated to a linear dependence of the 
temperature gradient on the diffuse fraction when Eqs. (A5) 
and (A6) are finally substituted into Eq. (3). 

The integral equations of the equivalent view factors 
over the surfaces of the beam are 


F ± (f m ;x i ) dxj = dxj 


(frn ; Xi , Vi , yj , Zj) dy, dyj dzj , 


(A7) 


F^{fm\Xi)dxi = dxi T\\(f m \Xi,yi,yj,zf)d.yidxjdyj , (A8) 


where the limits of the integrals depend on the surface m 
from which the phonons were originally radiated. 

Since a temperature gradient is only present in the prin- 
cipal direction of heat flow (x-axis), all surface elements at 
the plane of constant x are at the same temperature and radi- 
ate the same flux. Hence, the equivalent view factor for the 
direct (non-scattered) heat flux F D (f m ;xj) from the surface to 
the reference plane can be calculated analytically from Eq. 
(Al ) by considering the view factor of a infinitesimal strip of 
width w to the reference plane at x = 0. 32 
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